Relations Between Timing Features and Colors in Accreting 

Millisecond pulsars 



Steve van Straaten 1 , Michiel van der Klis 1 , Rudy Wijnands 

straatenOscience . uva.nl 



ABSTRACT 

We have studied the aperiodic X-ray timing and color behavior of the accreting millisecond 
pulsars SAX J1808.4-3658, XTE J1751-305, XTE J0929-314, and XTE J1814-338 using large 
data sets obtained with the Rossi X-ray Timing Explorer. We find that the accreting millisecond 
pulsars have very similar timing properties to the atoll sources and the low luminosity bursters. 
Based on the correlation of timing and color behavior SAX J1808. 4-3658 can be classified as an 
atoll source, and XTE J0929-314 and XTE J1814-338 are consistent with being atoll sources, 
but the color behavior of XTE J1751-305 is different. Unlike in other atoll sources the hard color 
in the extreme island state of XTE J 175 1-305 is strongly correlated with soft color and intensity, 
and is not anti-correlated with any of the characteristic frequencies. We found previously, that the 
frequencies of the variability components of atoll sources follow a universal scheme of correlations. 
The frequency correlations of the accreting millisecond pulsars are similar, but in the case of SAX 
J1808. 4-3658 and XTE J0929-314 shifted compared to those of the other atoll sources in a way 
that is most easily described as a shift in upper and lower kilohertz QPO frequency by a factor 
close to 1.5. Although, we note that the shift in lower kilohertz QPO frequency is based on only 
one observation for SAX J1808.4-3658. XTE J1751-305 and XTE J1814-338, as well as the low 
luminosity bursters show no or maybe small shifts. 
1. INTRODUCTION 

Many low mass X-ray binaries containing a neu- 
tron star show kilohertz QPOs, quasi periodic os- 
cillations with frequencies ranging from a few hun- 
dred to more than a 1200 Hz. The low frequency 
(< 200 Hz) part of the power spectrum taken from 
their X-ray lightcurves is usually dominated by a 
broad band- limited noise component. In addition 
to the band-limited noise several Lorentzian com- 
ponents below 200 Hz are present (see e.g. van 
Straaten et al. (2002)). These Lorentzians be- 
come broader as their characteristic frequency de- 
creases (Psaltis, Belloni, & van der Klis 1999a; 
van Straaten et al. 2002). Therefore, several fea- 
tures can, only at high frequencies, be classified 
as a QPO (traditionally defined by the condition 
FWHM < centroid frequency /2, see van der Klis 
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1995). In addition to these timing features in the 
persistent emission, 13 sources show millisecond 
oscillations (270-619 Hz) during thermonuclear X- 
ray bursts (see Muno 2004, for a review). Al- 
though these burst oscillations show a slight drift, 
they are thought to be connected with the spin of 
the neutron star and were the first direct indica- 
tion that these sources harbour neutron stars with 
millisecond spins. 

Based on the correlated behavior of their timing 
properties and the behavior of the sources in color- 
color diagrams, that are used to study their X-ray 
spectral properties, many of the neutron-star low- 
mass X-ray binaries are traditionally classified as 
either "atoll" or "Z" sources (Hasinger & van der 
Klis 1989). Recently, Muno et al. (2002) and in- 
dependently Gierlihski & Done (2002a) used the 
large Rossi X-ray Timing Explorer data sets now 
available for several of the atoll sources, to show 
that some of the atoll sources look much more like 
Z sources when plotted in color-color diagrams 



than was previously thought. Although the ap- 
pearances in the color-color diagram of the two 
classes now is more similar, the classification still 
holds, as clear differences between the two classes 
in both timing and color behavior remain (Muno 
ct al. 2002; Barret & Olive 2002; Olive et al. 2003; 
van Straaten et al. 2003; Reig et al. 2004). Note, 
that next to the differences in timing and color 
behavior there is also a difference in luminosity 
between the two classes; where the Z sources are 
all very luminous (<~ 3 x lO^ergs^ 1 Ford et al. 
2000), the atoll sources occur over a wide range of 
luminosities (10 36 ~ 38 ergs^ 1 Ford et al. 2000). 

The sources IE 1724-3045, GS 1826-24, and 
SLX 1735-269, are all low- luminosity neutron 
stars that show thermonuclear X-ray bursts. The 
power spectra of these "low-luminosity bursters" 
resemble those of the atoll sources in their hard 
states (Olive et al. 1998; Barret et al. 2000; Bel- 
loni, Psaltis, & van der Klis 2002); they show a 
broad band-limited noise component and several 
broad Lorentzian components. However, where 
the characteristic frequencies of the low-frequency 
Lorentzians were found to be similar to those in 
the atoll sources, the characteristic frequency of 
the highest frequency Lorentzian seemed to be 
higher (Belloni et al. 2002). 

In 1998 the first accreting millisecond pulsar 
was found with the Rossi X-ray Timing Explorer 
(RXTE) during an outburst of SAX J1808.4-3658 
(Wijnands & van der Klis 1998a). During this 
outburst no thermonuclear X-ray bursts were ob- 
served and the source did not show kilohertz QPOs 
above 200 Hz (Wijnands & van der Klis 1998b). 
The low frequency power spectrum was very sim- 
ilar to that of the atoll sources (Wijnands & van 
der Klis 1998b) and both the atoll sources and 
SAX J1808.4-3658 followed one relation when 
the characteristic frequency of a low-frequency 
Lorentzian was plotted versus that of the band- 
limited noise (Wijnands & van der Klis 1999). 
During a new outburst of SAX J1808.4-3658 in 
October 2002, kilohertz QPOs (up to 700 Hz) and 
burst oscillations were observed in an accreting 
millisecond pulsar for the first time (Wijnands 
et al. 2003; Chakrabarty et al. 2003). Recently, 
four additional accreting millisecond pulsars were 
discovered; XTE J1751-305, XTE J0929-314, 
XTE J1807-294, and XTE J1814-338 (see Wij- 
nands 2004, for a observational review). Of these 



XTE J1807-294 showed kiloherz QPOs (Mark- 
wardt, private communication) and XTE J1814- 
338 showed burst oscillations (Strohmayer et al. 
2003). 

Here we present an extensive analysis of the 
aperiodic-timing and color behavior of the accret- 
ing millisecond pulsars SAX J1808.4-3658, XTE 
J1751-305, XTE J0929-314, and XTE J1814-338. 
We study the correlations between the characteris- 
tic frequencies of their various timing features, and 
compare these with those of four well-studied atoll 
sources and three low- luminosity bursters. We 
find, that the accreting millisecond pulsars have 
very similar timing properties as the atoll sources 
and the low luminosity bursters, although dis- 
crepancies occur for SAX J1808.4-3658 and XTE 
J0929-314. We also compare the behavior of the 
accreting millisecond pulsars in the color-color dia- 
grams to that of the atoll sources and find that one 
pulsar, XTE J175 1-305, exhibits unusual behavior 
while the others act like ordinary atoll sources. In 
the next section we describe the observations and 
data analysis. We suggest that readers not inter- 
ested in the technical details of observations and 
analysis skip immediately to the results section. 

2. OBSERVATIONS AND DATA ANAL- 
YSIS 

In this work we study four accreting millisec- 
ond pulsars: SAX J1808.4-3658, XTE J1751-305, 
XTE J0929-314, and XTE J1814-338. We also 
re-analyzed data of Belloni et al. (2002) for the 
low luminosity bursters IE 1724-3045, GS 1826- 
24, and SLX 1735-269 (see §5.1). For the ac- 
creting millisecond pulsars we analyzed data from 
RXTE's proportional counter array (PCA; for in- 
strument information see Zhang et al. 1993); for 
SAX J1808.4-3658 we used data from March 11 to 
May 6, 1998 (the 1998 outburst) and from Octo- 
ber 15 to November 26, 2002 (the 2002 outburst), 
for XTE J 175 1-305 from April 4 to April 30, 2002, 
for XTE J0929-314 from May 2 to May 21, 2002, 
and for XTE J1814-338 from June 8 to August 
4, 2003. The data were divided into observations, 
identified by a unique RXTE observation ID, that 
cover one to several satellite orbits with ~ 3000 s 
of useful data each. We excluded data for which 
the angle of the source above the Earth limb is less 
than 10 degrees or for which the pointing offset is 
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greater then 0.02 degrees. We also excluded type I 
X-ray bursts. We found 4 bursts in the 2002 out- 
burst of SAX J1808.4-3658 (for results on these 
bursts see Chakrabarty et al. 2003) for which we 
excluded 16 s before and 384 s after the onset of 
each burst, 1 burst in XTE J1751-305 for which 
we excluded 16 s before and 592 s after the onset 
of the burst (note that this burst occurred during 
a weak reflare from April 27-30 for which it is un- 
clear whether it originated from XTE J 175 1-305, 
see Markwardt et al. 2002, the burst itself proba- 
bly originated from GRS 1747-312, which is how- 
ever ruled out as the source of the reflare, see in't 
Zand et al. 2003) and 22 bursts in XTE J1814-338 
(for results on these bursts see Strohmayer et al. 
2003) for which we excluded 16 s before and 240 s 
after the onset of each burst, the amount of time 
excluded depending on the length of the bursts. 
No bursts were found in XTE J0929-314. 

2.1. Color- Color Diagrams 

We used the 16-s time-resolution Standard 2 
mode to calculate the colors. For each of the 
PCA's five PCU detectors separately we calcu- 
lated a hard color, defined as the count rate in the 
energy band 9.7—16.0 keV divided by that in the 
6.0—9.7 keV band, and a soft color, defined simi- 
larly as 3.5-6.0 kcV/2.0-3.5 keV. For each detec- 
tor we also calculated the intensity, the count rate 
in the 2.0—16.0 keV band. To obtain the count 
rates in these exact energy ranges we interpolated 
linearly between count rates in the PCU channels 
based on the channel to energy calibration. We 
subtracted the background contribution in each 
band using 'pcabackest' version 3.0, with as input 
the "Mission-long Faint Model File" if the count 
rate in the full energy band becomes less than 40 
counts s^ 1 during an observation (faint data), and 
otherwise the "Mission-long Bright Model File" 
(bright data). We averaged up the background- 
subtracted 16-s count rates to obtain count rates 
per observation, and we calculated the colors sep- 
arately for each PCU. An observation spanned be- 
tween 1 and 30 ks, and within an observation the 
colors generally varied by 3 to 9 times the 1 a stan- 
dard deviation of the 16-s colors. For faint data 
we filtered out data taken within thirty minutes of 
the peak of the South Atlantic Anomaly and data 
with high electron contamination. 

The voltage settings of the PCUs on board 



RXTE were changed on three occasions defining 
the first four RXTE gain epochs. The fifth RXTE 
gain epoch started when a micro-meteorite created 
a small hole in the front window of PCU0. In order 
to correct for the changes in effective detector area 
between the different RXTE gain epochs and for 
the gain drifts within these epochs, as well as for 
the differences in effective area between the PCUs 
themselves we calculated the colors for each PCU 
of the Crab nebula, which can supposed to be con- 
stant in its colors, in the same manner. We then 
averaged the 16-s Crab colors and intensity per 
PCU for each day. For each PCU we divided the 
obtained color and intensity points per observation 
by the corresponding Crab values that are closest 
in time but within the same gain epoch. We then 
averaged the colors and intensity over all PCUs. 
Note that starting May 12, 2000, the propane layer 
on PCU0, which functions as an anti-coincidence 
shield for charged particles, was lost, defining the 
beginning of gain epoch 5. Therefore we applied 
the additional filter recommended by the RXTE 
science team to PCU0 data from epoch 5. 

2.2. Timing 

For the Fourier timing analysis we used the 
122— [is time-resolution Event modes and the 
0.95— /us time-resolution Good Xenon modes. We 
used data from all energy channels. We con- 
structed power spectra per observation (see above) 
using data segments of 256 s and 1/8192 s time 
bins such that the lowest available frequency is 
1/256 Hz and the Nyquist frequency 4096 Hz; the 
normalization of Leahy et al. (1983) was used. We 
subtracted a Poisson noise level using the method 
of Klein- Wolt et al. (2004), which is built on the 
analytical function from Zhang et al. (1995). The 
resulting power spectra were then converted to 
squared fractional rms. Within the observations 
the power spectra did not change within the er- 
rors. When possible several observations, adjacent 
in time, and for which the power spectra remain 
the same, were added up to improve the statistics. 
These are called groups. Occasionally, several ob- 
servations, not adjacent in time, for which the 
power spectra are the same, were added up into 
groups. The groups for which this was done will 
be noted in the text (see §3). Before fitting, all 
frequency bins (generally about 20) containing the 
pulse spike were removed from the power spectra. 
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As a fit function we used the multi-Lorentzian 
function; a sum of Lorentzian components (see e.g. 
Belloni et al. 2002; van Straaten et al. 2002). Note, 
that we no longer constrained the Lorentzian (van 
Straaten et al. 2002, 2003) that fits the band- 
limited noise component to be zero-centered, but 
instead let all Lorentzian parameters float. This 
was motivated by the work of Pottschmidt et al. 
(2003), who found that for the band-limited noise 
component in the black hole candidate Cyg X— 1, a 
non zero-centered Lorentzian provided better fits 
(usually we find a Q significantly larger than zero 
but < 0.3 for this component, see below). We 
generally only included those Lorentzians in the 
fit whose power can be measured to an accuracy 
of better than 33%. The few exceptions where 
we deviated from this will be noted in the text 
(§3). We plot the power spectra in the power times 
frequency representation (yV v ), where the power 
spectral density P„ is multiplied by its Fourier 
frequency v. For a multi-Lorentzian fit function 
this representation helps to visualize a character- 
istic frequency, ^ max , the frequency where each 
Lorentzian component contributes most of its vari- 
ance per logarithmic frequency interval, as in vY v 
the Lorentzian's maximum occurs at ^ max ("max = 
\Jvq + A 2 , where f is the centroid and A the 
HWHM of the Lorentzian; Belloni ct al. (2002)). 
We represent the Lorentzian relative width by Q 
defined as ^o/2A. 

3. RESULTS 

We find that three to seven Lorentzian com- 
ponents arc needed to fit the power spectra of 
the millisecond pulsars. All Lorentzian compo- 
nents can be identified in the low mass X-ray 
binary scheme of Belloni ct al. (2002) and van 
Straaten et al. (2003). In Figure 1 we illustrate the 
main multi-Lorentzian components of this scheme. 
Note, that some updates to the scheme based on 
results of this paper are made (see §5.2). The fit 
parameters are listed in Tables 1 (^ ma x), 2 (Q), 
and 3 (fractional rms amplitude) . Of the fits for 
the accreting millisecond pulsars listed in Tables 
1, 2, and 3 47% has a x 2 /dof below 1.1, 27% has a 
X 2 /dof between 1.1 and 1.3, 23% between 1.3 and 
1.7, and 3% (group 4) has X 2 /dof = 2.2. The de- 
grees of freedom range between 91 and 142. Note, 
that the x 2 /dof values of the multi-Lorentzian fit 
function are quite high. However, they are com- 



parable to previous results of this fit function (van 
Straaten et al. 2002, 2003), and other fit func- 
tions tried in previous studies (see Belloni et al. 
2002; van Straaten et al. 2002, for a comparison 
between the multi-Lorentzian fit function and the 
broken powerlaw fit function ) do not give better 
results. Note that this is not a problem as here 
we are only measuring the relevant characteristic 
frequencies and not trying to find the true model 
for the power spectra. The fit groups contain data 
from several RXTE observations. The observation 
IDs of each fit group are listed in Table 4. 

3.1. SAX J1808.4 3658: the 2002 Out- 
burst 

In this paragraph we describe the color and tim- 
ing behavior of SAX J1808.4-3658 during the 2002 
outburst as a function of time with the help of Fig- 
ures 2 and 3. In Figure 2 we show six representa- 
tive power spectra for the 2002 outburst. In Figure 
3 we plot hard color, soft color, intensity (see also 
Wijnands et al. 2003) and the frequency of one of 
two timing components (depending on epoch) ver- 
sus time. We use the terminology of Belloni et al. 
(2002) and van Straaten et al. (2003) in identifying 
the timing features (each Lorentzian component is 
indicated with an L, for example the upper kilo- 
hertz QPO is called L M ). In the first group (group 
1) a broad L& at ~ 10 Hz, Lh at ~ 50 Hz, L^Hz 
at ~ 200 Hz and L u at ~ 600 Hz are observed (see 
Fig. 2). As the intensity increases to a maximum, 
both the hard color and the soft color drop to a 
minimum (Fig. 3). The power spectra (groups 
2 and 3) now show L u at a higher characteristic 
frequency than in group 1. LhHz is also present 
but at a rather high frequency (~ 300 Hz). The 
low frequency part of the power spectrum is com- 
plex; L/j is at a higher frequency than in group 
1, and below L^ another low frequency QPO ap- 
pears that can be identified with Llf (see §5.2). 
Lf, changes from a broad band-limited noise com- 
ponent into a QPO and Lfc2 appears at a lower 
characteristic frequency (see §5.2). Group 3 is the 
only group where the lower kilohertz QPO, Lg, 
can be detected (see also Wijnands et al. 2003). 
In the fit we fixed the Q value to the value that 
was found by Wijnands et al. (2003). This was 
done because with the rebinning used in our fit 
we only had a few points in this peak. After the 
maximum, the intensity starts to decrease while 
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the hard color and the soft color increase. The 
power spectrum of group 4 is similar to that of 
group 1, showing L&, Lh, LhHz, and L u , at some- 
what lower characteristic frequencies. While the 
intensity continues to decrease, the power spec- 
trum (group 5) still shows these components with 
Vb, Vhi an d v u even lower. The characteristic fre- 
quency of L u is now below the spin frequency of 
401 Hz. While the intensity keeps decreasing and 
the soft color starts decreasing, the power spec- 
trum remains about the same for groups 5—8. 
After that there is a slight increase in intensity 
accompanied by a decrease in hard color and an 
increase of Vb, Vh, and v u (group 9). Then, while 
intensity and soft color decrease again, the power 
spectra of groups 10—15 are very similar to that 
of group 5. It is not influenced by the break in the 
intensity and soft color versus time curves that oc- 
curs at MJD 52575. In group 15, LhHz is no longer 
observed with an upper limit consistent with this 
component still being present at similar strength 
(see Table 3). After that no power spectral com- 
ponents are detected until MJD 52578 due to bad 
statistics (see Wijnands et al. 2003). After MJD 
52578 the behavior of the source has changed com- 
pletely; the intensity fluctuates on time scales of 
days, and the power spectra show a strong QPO 
around 1 Hz (see Fig. 2). This behavior is similar 
to that observed during what appeared to be the 
end stage of an outburst in 2000 (Wijnands et al. 
2001; Wijnands 2004). The 1 Hz QPOs could not 
be fitted well with Lorentzians. Instead, several 
Gaussians were required. In addition to the 1 Hz 
QPO sometimes another QPO is present at around 
30 Hz. As an illustration the centroid frequency 
of the main Gaussian used to fit the 1 Hz QPO 
is plotted in the bottom panel of Figure 3. The 
full behavior of the 1 Hz QPO is beyond the scope 
of this paper and a detailed analysis of this phe- 
nomenon is in progress. 

During the 2002 outburst Wijnands et al. 
(2003) detected a narrow QPO near 410 Hz. To 
search for these 410 Hz QPOs we created power 
spectra similar to those of §2.2 but now in the 
3-60 keV energy range. This is the same en- 
ergy range for which Wijnands et al. (2003) found 
these QPOs. We removed the Poisson level and 
the pulse spike and fitted the power spectra with 
the multi-Lorentzian fit function (see also §2.2). 
To be able to detect the narrow 410 Hz QPO we 



had to use a finer rebinning than in §2.2. We 
sometimes had to add several groups, adjacent in 
time, and for which the power spectra remained 
the same, to detect the 410 Hz QPO. We find 
the 410 Hz QPO in group 9 and in the combined 
groups 5-6, and 10-11. In Figure 4 we plot the 
combined power spectrum of groups 5-6, 9 and 
10-11. The characteristic frequency of the 410 Hz 
QPO, the frequency difference between the ~ 401 
Hz pulse spike and the 410 Hz QPO, and Lh are 
listed in Table 5. We do not detect 410 Hz QPOs 
during the 1998 outburst. 

3.2. SAX J1808.4 3658: the 1998 Out- 
burst 

In this paragraph we describe the color and 
timing behavior of SAX J1808.4-3658 during the 
1998 outburst as a function of time with the help 
of Figures 5 and 6. The 1998 outburst of SAX 
J1808. 4-3658 has been studied extensively pre- 
viously; the broad-band timing properties where 
studied by Wijnands & van der Klis (1998b), and 
the X-ray energy spectrum was studied by Gil- 
fanov et al. (1998), Heindl & Smith (1998), and 
Gicrlihski et al. (2002). In group 16 the power 
spectrum shows Lb at ~ 2 Hz, Lh at ~ 10 Hz, LhHz 
at - 190 Hz and L u at - 270 Hz (see Fig. 5). The 
characteristic frequencies are all lower than during 
the 2002 outburst. Note, that the power of LhHz 
and L u could only be measured to an accuracy of 
better than 50% (and not 33%, see §2.2). However, 
omitting either one of these components leads to 
one broad Lorentzian fitting both LhHz and L u . 
For the second observation (at MJD 50916) there 
is not enough data left to create a power spectrum 
after the screening criteria have been applied (see 
§2). After a gap of four days, the intensity and 
soft color are lower and the hard color is higher. 
Then, as the intensity and soft color decrease, the 
power spectrum and hard color remain the same 
during groups 17—20. The power spectrum (see 
group 19 in Fig. 5) shows Lb, Lh, and L u at the 
lowest characteristic frequencies observed for SAX 
J1808. 4-3658. It also shows Lg ow at about 30 Hz. 
In group 21 the intensity, soft and hard color are 
lower and Vb, Vh, vi ow , and v u are higher (see Fig. 
5). At MJD 50929, a change of slope in the inten- 
sity and soft color curves occurs. This change of 
slope was also observed during the 2002 outburst 
(see Fig. 6). Group 22 was fitted with only two 
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Lorentzians, L;, and a broad Lorentzian fitting the 
higher frequencies. Due to lack of statistics L/,, 
Leow, and L u could not be identified. 

3.3. XTE J1751 305 

In this paragraph we describe the color and tim- 
ing behavior of XTE J 175 1-305 during its 2002 
outburst as a function of time with the help of Fig- 
ures 7 and 8. In Figure 7 we show three representa- 
tive power spectra. In Figure 8 we plot hard color, 
soft color, intensity and versus time. Note that 
in our analysis we subtract the background with 
RXTE background models (see §2.1). In addition 
to this there is also a significant background com- 
ponent for XTE J 175 1-305 due to nearby sources 
and Galactic diffuse emission (Markwardt et al. 
2002). We draw the baseline (the dashed line in 
the insert of the intensity curve in Fig. 8) through 
the same observations as was done by Markwardt 
et al. (2002) and omit the soft and hard color 
points that are obtained when the source is ac- 
tually undetected. 

In the first observation (group 23) XTE 
J1751-305 shows L& at ~ 1 Hz, L/j at ~ 6 Hz, 
Lfoto at ~ 30 Hz, and L u at ~ 280 Hz (see Fig. 
7). Then with intensity, hard color and soft color 
decrease (see Fig. 8), group 24 shows the same 
power spectral components as group 23, but at 
lower characteristic frequencies. One day later 
the characteristic frequencies are even lower. Then 
the power spectrum remains the same for 6 days. 
To further improve the power spectral statistics 
we combined these 6 days into group 25. At 
MJD 52575, the intensity, hard and soft color 
curves show a change in slope similar to that in 
the intensity and soft color curves of the 1998 and 
2002 outbursts of SAX J1808.4-365 (see §3.1 and 
3.2). After the observations that span up group 
25, the power spectrum still shows power but the 
statistics are not sufficient to identify the power 
spectral components. Between MJD 52380 and 
MJD 52391 the source is not detected, then a 
small reflare is observed (see the insert of the in- 
tensity curve in Fig. 8), for which it is uncertain 
whether it is from XTE J 175 1-305 (Markwardt 
et al. 2002). The statistics are too low to detect 
any power spectral components during the reflare. 



3.4. XTE J0929 314 

In this paragraph we describe the color and tim- 
ing behavior of XTE J0929-314 during its 2002 
outburst with the help of Figures 9 and 10. In 
Figure 9 we show three representative power spec- 
tra. In Figure 10 we plot hard color, soft color, 
intensity and if, versus time. For the first obser- 
vation (at MJD 52396) there are not enough data 
left to create a power spectrum after the screening 
criteria have been applied (see §2). Then, after 
a gap of seven days, the intensity and soft color 
are lower and the hard color is higher (see Fig. 
10). As the intensity and soft color keep declining 
the hard color becomes approximately constant. 
In the next seven days (MJD 52403-52410) the 
power spectrum shows at low frequencies L&, L/,, 
and Ijlf/2- is a narrow QPO {Q « 3) with 

a characteristic frequency that varies between 0.6 
and 1.3 Hz. At high frequencies the power spec- 
trum shows one broad Lorentzian that probably 
(see below) fits both Le ow and L„. To improve 
the statistics at high frequencies we add up the 
observations with 0.64 < Vlf/2 < 0.86 into group 
26 and the observations with 0.99 < vlf/2 < 1-33 
into group 27. The power spectrum of group 26 
shows Lb at ~ 0.25 Hz, L/j at ~ 2 Hz, and L LF / 2 at 
~ 0.7 Hz (Fig. 9). At high frequencies Le ow and 
L u could now be identified. Group 27 is simi- 
lar to group 26, except that, as mentioned above 
at high frequencies the fit could not distinguish 
between Li ow and L u (Fig. 9). After that, from 
MJD 52411-52428, the power spectra still show 
power but the statistics are not sufficient to iden- 
tify the power spectral components. If we add up 
all observations from MJD 52411-52428 we obtain 
group 28. This group is similar to group 27, and 
just as in group 27 the fit does not distinguish be- 
tween h low and L u (Fig. 9). After MJD 52428 the 
statistics become insufficient and the power spec- 
tra no longer show significant power. 

3.5. XTE J1814 -338 

In this paragraph we describe the color and tim- 
ing behavior of XTE J1814-338 during its 2003 
outburst as a function of time with the help of 
Figures 11 and 12. In Figure 11 we show three 
representative power spectra. In Figure 12 we plot 
hard color, soft color, intensity and Vb versus time. 
From MJD 52798 to MJD 52825 the power spectra 
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at low frequencies show Lb at ~ 0.5 Hz and at 
~ 3 Hz. The characteristic frequencies vary sys- 
tematically (see Fig. 12), i/j between 0.25 and 0.61 
Hz, and Vh, which correlates with i>b, between 1.4 
and 3.7 Hz. At high frequencies the power spectra 
are fitted with one broad Lorentzian that proba- 
bly (see below) fits both Li ow and L u . To improve 
the statistics at high frequencies we add up the 
observations with 0.31 < Vb < 0.39 into group 29 
and the observations with 0.41 < v\> < 0.61 into 
group 30. Group 29 shows Lb at ~ 0.38 Hz, 
L/j at ~ 2.2 Hz, Iji ow Hz at ~ 21 Hz and L„ at 
~ 190 Hz (see Fig. 11). It also shows the nar- 
row L iF / 2 at ~ 0.9 Hz. Group 30 shows, except 
for L LF / 2 > the same components as group 29 (see 
Fig. 11). In group 30 all characteristic frequencies 
are higher; Vb w 0.5 Hz, Vh w 3 Hz, vg w 50 Hz, 
and v u w 220 Hz. After MJD 52825 the power 
spectra still show power between MJD 52826 and 
MJD 52829, but the statistics are not sufficient 
to identify the power spectral components. Af- 
ter MJD 52830 the statistics become insufficient 
and the power spectra no longer show significant 
power. 

4. Timing and Colors 

Now that we have described and identified all 
power spectral components in the accreting pul- 
sars SAX J1808.4-365, XTE J1751-305, XTE 
J0929-314, and XTE J1814-338, we can link the 
timing properties to the position in the color-color 
diagram. We can also compare the timing prop- 
erties and X-ray spectral properties to the classic 
island and banana states described by Hasingcr 
& van der Klis (1989) and compare with other 
sources. An addition to the classic states described 
by Hasinger & van der Klis (1989) is the extreme 
island state. At the lowest count rates, an exten- 
sion of the island state is sometimes seen with a 
harder spectrum and stronger band-limited noise 
than the "canonical" island state. The term "ex- 
treme island state" has been used to designate 
such a state (Prins & van der Klis 1997; Reig et al. 
2000). In Reig et al. (2000) this extreme island 
state was detached from the "normal" island state. 
For 4U 1608-52 (van Straaten et al. 2003) and 
now here for SAX J1808.4-3658 the change from 
extreme to " normal" island state is more gradual, 
so the question arises where to put the boundary. 
Here we have used the presence of Lg ow . In the ex- 



treme island state the power spectrum shows he w 
and in the island state Lg ow is absent. 

4.1. SAX J1808.4 3658 

In Figure 13 we plot the color-color and soft 
color vs. intensity diagram of both the 2002 
and 1998 outburst of SAX J1808.4-365. For the 
2002 outburst we only plot observations from be- 
fore MJD 52578, at which time the behavior of 
the source changed completely (see §3.1). At 
the start of the 2002 outburst (group 1) the 
power spectrum is that of an island state; it shows 
strong broadband noise with low ^ max and no 
very low frequency noise (VLFN). Group 2 and 
3 show several power spectral characteristics of 
the lower left banana, the characteristic frequen- 
cies are higher, L£ just appears, and the low fre- 
quency part of the power spectrum becomes com- 
plex. Only the VLFN that is typical for the ba- 
nana state has not yet appeared. In group 4 the 
source is back in the island state. After that the 
source remains in an island state. The charac- 
teristic frequencies decrease until group 5, and 
after that only small changes in the character- 
istic frequencies occur. During the whole 2002 
outburst SAX J1808. 4-365 did not go into an ex- 
treme island state. In an extreme island state the 
power spectra show Lb, L/,, Lg ow and L u all at 
very low characteristic frequencies (van Straaten 
et al. 2003). In 4U 1608-52 (van Straaten et al. 
2003) the boundary between island and extreme 
island state lies at a Vb of about 2-3 Hz. In the 
2002 outburst of SAX J1808.4-3658 Uow was not 
observed and Vb was always above 2 Hz. During 
the 1998 outburst the characteristic frequencies of 
all components were lower than during the 2002 
outburst. In the first group (group 16) of the 
1998 outburst SAX J1808.4-365 was in the island 
state. After that the source was in an extreme 
island state showing Li ow and Vb < 1. 

The determination of the atoll source states 
based on the power spectral results is in agreement 
with the positions in the color-color diagram. SAX 
J1808. 4-365 behaves similarly to the atoll sources 
4U 1705-44, Aql X-l, and 4U 1608-52 (Muno 
et al. 2002; Gierlihski & Done 2002a; Barret & 
Olive 2002; Olive et al. 2003; van Straaten et al. 
2003; Reig et al. 2004) which in the extreme island 
state show horizontally extended branches in the 
color-color diagram. At the beginning of the 2002 
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outburst SAX J1808. 4-365 showed a state transi- 
tion from the island state (group 1) to groups 
2 and 3 and back to the island state (group 
4). This transition is seen in the color-color di- 
agram as a lowering of the hard and soft color 
and a slight increase in intensity. At the lowest 
hard color (groups 2 and 3) the source has almost 
reached the lower left banana (see above). It seems 
unlikely that the actual banana branch is reached 
as the gap between group 2 and 3 is only 3000 s. 
Then the hard color increases again and the source 
goes back to an island state. As noted previously 
(Mendez et al. 1999; van Straaten et al. 2003), 
hard color appears to be the parameter best cor- 
related to the power spectral characteristics in the 
island states, and this is confirmed here. During 
the whole transition, and also in the other groups, 
the hard color is anti-correlated with the charac- 
teristic frequencies of the power spectral compo- 
nents (see Figure 14). This behavior is typical 
for atoll sources (4U 0614+09, van Straaten et al. 
2000; 4U 1728-34, Di Salvo et al. 2001; 4U 1705- 
44, Olive et al. 2003; 4U 1608-52, van Straaten 
et al. 2003). Note that compared to other atoll 
sources the relative change in hard color during the 
transition is only subtle. For SAX J1808.4-3658 
the relative change in hard color is about 5% were 
it is > 30% for 4U 1728-34 (Di Salvo et al. 2001), 
4U 1705-44 (Olive et al. 2003), and 4U 1608-52 
(van Straaten et al. 2003). The relative change in 
soft color is about 10%, comparable to the rela- 
tive soft color changes for 4U 1728-34, 4U 1705- 
44, and 4U 1608-52. The behavior after the state 
transition can be described in the scheme pro- 
posed for the atoll sources by van Straaten et al. 
(2003) (see their figure 12). The source remained 
in an island state in which the characteristic fre- 
quencies arc low and hardly change after group 
5. Within the island state the soft color changes 
in correlation with intensity. As the intensity de- 
cays the soft color decreases and a horizontally 
extended island branch appears. Note, that the 
constant frequencies make this branch very differ- 
ent from the Z source horizontal branch. In other 
sources horizontally extended branches where only 
observed in the extreme island state (van Straaten 
et al. 2003), here we observe the same phenomenon 
but now in an island state. However, note that 
the characteristic frequencies in this island state 
are low and only the absence of Ij£ OW distinguish 



it from an extreme island state (see above) The 
small changes in characteristic frequency that do 
occur are accompanied by anti-correlated changes 
in hard color (see Fig. 14) and show no relation 
to soft color or intensity. In the 1998 outburst the 
same phenomenon occurs; in groups 17—20 the 
power spectra remain in one extreme island state 
and as the intensity decreases so does the soft color 
and again a horizontally extended extreme island 
branch appears. In groups 21—22 the character- 
istic frequencies are higher and another horizon- 
tally extended extreme island branch is formed at 
a lower hard color. Thus in total three horizontally 
extended island branches (5-15, 21-22, and 17-20) 
appear above each other in the left panel of Fig- 
ure 13, and the frequencies on each of them anti- 
correlate with hard color. This anti-correlation 
with hard color also holds between the 1998 and 
the 2002 outbursts. For a large part (from group 
5 of the 2002 outburst) both outbursts cover the 
same intensity range (see Fig 13). At similar in- 
tensities the characteristic frequencies in the 1998 
outburst are much lower, while the hard color is 
higher. 

4.2. XTE J1751 305 

In Figure 15 we plot the color-color and soft 
color vs. intensity diagram of XTE J1751-305. 
The power spectra of XTE J1751-305 show Lb, L/j, 
hi ow and L„ all at low characteristic frequencies. 
This is typical for an extreme island state. The 
source stays in an extreme island state during the 
whole outburst while the characteristic frequencies 
change between v\> = 1.2 and 0.4 Hz. Although 
the power spectra are typical for those of the ex- 
treme island of the atoll sources, the behavior in 
the color-color diagram is not. During group 25 
the characteristic frequencies remain the same. In 
the atoll sources (and in SAX J1808. 4-365, see 
above) the hard color would then remain the same 
while the intensity and soft color change in a cor- 
related fashion. This would lead to an extended 
horizontal branch in the color-color diagram. Here 
the hard color is not constant but is correlated 
with intensity in a similar way as the soft color. 
This correlation between hard color, soft color and 
intensity exists during the whole outburst. Also, 
unlike in the atoll sources and SAX J1808. 4-3658, 
no anti-correlation occurs between hard color and 
the characteristic frequencies of the power spectral 
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components (see Fig. 14). In Figure 15 we also 
plot the colors of the reflare and those of the back- 
ground due to nearby sources and Galactic diffuse 
emission (see §3.3). The points of the reflare do 
not lie on the the curve of the main outburst, but 
on an extension of the background. This suggests 
that the reflare is caused by one of the background 
sources and not by XTE J1751-305. 

4.3. XTE J0929 314 

In Figure 16 we plot the color-color and soft 
color vs. intensity diagram of XTE J0929-314. 
The filled squares mark group 26, the open dia- 
monds group 27, and the filled stars group 28. The 
source remains in an extreme island state during 
the whole outburst, showing Lb, L/j, Le ow and L u 
all at low characteristic frequencies. The power 
spectra also show the narrow Llf/2- This com- 
ponent is also found in 4U 1608-52 (van Straaten 
ct al. 2003) and in XTE J1814-338 (see §4.4), also 
in the extreme island state. So, the timing be- 
havior is that of the extreme island state of the 
atoll sources. The behavior in the color-color dia- 
grams is also similar to that of an extreme island 
state. Within the extreme island state the soft 
color changes in correlation with intensity. As the 
intensity decreases the soft color decreases and a 
horizontally extended island branch appears (see 
Fig. 16). The hard color and the characteristic fre- 
quencies of the power spectral components are not 
correlated with intensity (and thus soft color). In 
the atoll sources the hard color is anti-correlated 
with the characteristic frequencies of the power 
spectral components. Here that anti-correlation 
is not evident, but the changes in hard color and 
characteristic frequencies are only small (see Fig. 
14). 

4.4. XTE J1814-338 

In Figure 17 we plot the color-color and soft 
color vs. intensity diagram of XTE J 1814-338. 
The filled squares mark group 29, and the open 
diamonds group 30. The source remains in an ex- 
treme island state during that part of the outburst 
for which power spectral fits could be made. It 
shows Lb, L/j, Le ow and L u all at low character- 
istic frequencies. Group 29 also shows the nar- 
row L LF / 2 - This component was also found in 
4U 1608-52 (van Straaten et al. 2003) and XTE 
J0929-314 (see §4.3), also in the extreme island 



state. The timing properties could only be mea- 
sured in a small range of the color-color diagrams 
(see Fig. 17). Within this range, the characteris- 
tic frequencies, colors and intensity only vary by a 
small amount, and therefore it is hard to compare 
this extreme island state of XTE J1814-338 with 
that of the atoll sources. Within the extreme is- 
land state of XTE J1814-338 the soft color (and 
intensity) does not change enough for the source 
to draw up a horizontal branch as is often seen 
in atoll sources in the extreme island state. The 
range in characteristic frequencies within the ex- 
treme island state is too small to observe the anti- 
correlation between hard color and the character- 
istic frequencies of the power spectral components 
that is seen in the atoll sources (see Fig. 14). 

5. DISCUSSION 

5.1. Correlations of Timing Features 

The frequencies of the variability components 
of the atoll sources, 4U 0614+09, 4U 1608-52, and 
4U 1728-34, follow a universal scheme of correla- 
tions when plotted versus v u . In Figure 18 we 
present such a plot for SAX J1808.4-3658, XTE 
J1751-305, XTE J0929-314, and XTE J1814-338 
(colored points). We also include (black points) 
results for the atoll sources 4U 0614+09, 4U 1728- 
34 (van Straaten et al. 2002), 4U 1608-52 (van 
Straaten et al. 2003), and Aql X-l (Reig et al. 
2004) . For Aql X-l we only include the results for 
the extreme island state for which v u could unam- 
biguously be identified. For 4U 1728-34 we also 
include a recent result by Migliari, van der Klis, 
& Fender (2003), where for the first time twin kilo- 
hertz QPOs were observed in the island state. As 
a third group of sources we have included the low 
luminosity bursters IE 1724-3045, GS 1826-24, 
and SLX 1735-269 (grey points) described in Bel- 
loni et al. (2002). Belloni et al. (2002) constructed 
power spectra only up to 512 Hz, which led to an 
ill-constrained L„. Here we re- analyze the timing 
properties of groups C-R of Belloni ct al. (2002) 
using the same method as described in §2.2 (Tables 
1, 2, and 3). With 0.0039-4096 Hz power spectra 
this leads to a much better constrained L M . 

SAX J1808. 4-3658 shows similar relations be- 
tween its characteristic frequencies as the atoll 
sources do. However, with respect to these 
sources, the relations of SAX J1808. 4-3658 are 
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shifted. XTE J0929-314 shows a similar shift, al- 
though this is only based on one power spectrum. 
XTE J1751-305 and XTE J1814-338 as well as the 
low luminosity bursters show relations consistent 
with those of the atoll sources. The shift occurs 
only between the characteristic frequencies of the 
low frequency components on the one hand and 
v u (and probably vg) on the other (see below): 
the mutual relations between the characteristic 
frequencies of the low frequency components are 
identical for all sources (see Figure 19). 

To determine the actual shift factors between 
the frequency relations of SAX J1808. 4-3658, 
XTE J0929-314, and those of the atoll sources 
we consider the Vi vs. v u relations, where Vi is 
either v^, Vh, or vg ow . We only use that part of 
the relations for which v u < 600 Hz, as the behav- 
ior of the low-frequency components above 600 Hz 
is complex (see also §5.4). We take a i/j vs. v u 
relation of one source (for example SAX J1808.4- 
3658) together with that of the atoll sources and 
perform a power law fit. As our points have errors 
in both coordinates, we use the FITEXY routine 
of Press Press & Tcukolsky (1992) that performs 
a straight-line fit to data with errors in both coor- 
dinates. We take the logarithm of both Vi and v u 
so that fitting a straight-line becomes equivalent 
to fitting a power law. Before fitting we multiply 
v u (or Vi) with a shift factor. We let this shift 
factor run between 0.1 and 3 with steps of 0.001. 
The fit with the minimal \ 2 then gives the shift 
factor. The errors in the shift factor use A% 2 = 
1.0 (corresponding to a 68% confidence level). In 
Table 6 we list the shift factors both in v u and 
v i for the different relations and sources. Next to 
the sources that show a clear shift in Figure 18 
(SAX J1808.4-3658 and XTE J0929-314) we also 
include the other millisecond pulsars as well as the 
low luminosity bursters. We also determine the 
shift factor for the single vg observation of SAX 
J1808.4-3658. We determine the shift factor if 
only vt is shifted, and the shift factor if both vg 
and v u arc shifted. For the latter we first multiply 
v u with 1.454, the weighted average of the v u shift 
factors for the Vb, Vh, and vg ow vs. v u relations of 
SAX J1808.4-3658 (see below). 

The shift factors in v u for SAX J1808.4-3658 
are between 1.4 and 1.5 when determined with 
the vi> and vt vs. v u relation. For the vg ow vs. 
v u relation the factor is higher, but this is based 



on fewer points which results in a larger error on 
the shift factor. From this point on we take the 
shift in v u to be the weighted average of the three 
which results in 1.454±0.009. After shifting v u 
with this factor, the additional shift needed for 
the single vg point is 1.359, similar to the shift 
factors found for v u . For XTE J0929-314 the v u 
shift factors determined with the Vb and Vh vs. 
v u relation are similar to those in SAX J1808.4- 
3658, but note that these results are based on only 
one point from XTE J0929-314. The v tow vs. v u 
relation, just like for SAX J1808.4-3658, gives a 
higher value but with large error. For some of the 
other sources small shifts might exist, but unlike 
for SAX J1808. 4-3658 the frequency vs. frequency 
relations for these sources have up to know only 
been observed over a small frequency range. The 
frequency relations should be measured over larger 
frequency ranges to confirm or reject the idea that 
small shifts occur. 

The simplest explanation for this shift is that 
some physical difference between these sources af- 
fects v u and vg. To make the relations of SAX 
J1808.4-3658 and XTE J0929 314 coincide with 
those of the atoll sources, v u and vg must be multi- 
plied with ~ 1.5. In Figure 20 we replot Figure 18, 
but now with the v u and v t of SAX J1808.4-3658 
and XTE J0929-314 multiplied by 1.454. The al- 
ternative, that vu and the characteristic frequen- 
cies of all the low frequency components, L;,, L/,, 
Lg ow , and the narrow QPOs (see §5.2 for a dis- 
cussion about the narrow QPOs) are affected is 
more complex. In this case to make the relations 
coincide the characteristic frequencies of all the 
low frequency components have to be multiplied 
with about 0.35, vg would have to be multiplied 
with about 0.8, and v^Hz would have to remain 
the same. It is simpler if v u and vg are the shifted 
frequencies. Note, that the resulting x 2 /dof of the 
power law fits (see Table 6) are independent of 
whether v u or the frequencies of the low frequency 
components are shifted. 

We tried to further test our idea that the shift 
occurs in v u and not in the frequencies of all the 
low frequency components, by plotting Q versus 
^max for both L u and L/j (see Figure 21). The 
Q u vs. v u plot showed a shift in v u similar to 
that in Figure 18, while the Qh vs. Vh plot did 
not show a clear shift, consistent with the sim- 
pler interpretation above. However, the scatter in 
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the Q versus i-Wx correlations is much larger than 
of the frequency-frequency correlations in Figure 
18. Neither plots of fractional rms versus y mM nor 
linking Figure 18 to the atoll source states led to 
extra clues. 

Might we have wrongly identified the high- 
est frequency Lorentzian component as L u while 
it really is L{t There are a number of argu- 
ments against this, all based on similarities with 
the atoll source class of which SAX J1808.4-3658 
seems part (see §4.1). First, Li was detected to- 
gether with the component identified as L u in SAX 
J1808.4-3658 by Wijnands et al. (2003) (our group 
3, see §3.1) when the source was almost in the 
lower left banana (see §4.1). In group 3 Lu is weak 
and L u is strong (see Table 3); this is common 
behavior for atoll sources when they are almost 
in the lower left banana (Mendez et al. 2001; van 
Straaten et al. 2002, 2003). Second, in the island 
and extreme island states of SAX J1808. 4-3658 
the component identified as L u is a strong broad 
component that is present over a broad frequency 
range while Lg is not. This is again common be- 
havior for the atoll sources, where Li was observed 
in an island state only once, and then as a narrow 
and weak QPO (Migliari et al. 2003). 

5.2. Identifying Timing Features 

The frequency relations of the accreting mil- 
lisecond pulsars and the low luminosity bursters of 
Figure 18, and especially the shift between the re- 
lations of SAX J1808. 4-3658 and the atoll sources, 
provide new clues about the identification of some 
power spectral components. We assume here that 
the shift in Figure 18 is a shift in v u ; the con- 
clusions drawn below do not change if we use the 
alternative, a shift in vi OWl Vh, an d Vb (see §5.1). 

The strong, broad highest-frequency Lorentzian 
in the extreme island state of the atoll sources 
could previously be identified as either L u or LhHz 
based on its frequency of about 200 Hz. The rms 
fractional amplitude vs. frequency relations of 
4U 0614+09, 4U 1608-52, and 4U 1728-34 gave 
a hint that this Lorentzian was L u (van Straaten 
et al. 2003). Now we observe the same shift in 
frequency for this Lorentzian as for the compo- 
nent which in other states is unambiguously L u 
(see Fig. 18), supporting the conclusion that it 
is indeed L u . This conclusion is further strength- 
ened by the results for the low luminosity bursters 



IE 1724-3045, GS 1826-24, and SLX 1735-269. 
These sources, which always have a power spec- 
trum similar to that in the extreme island state of 
the atoll sources, show a clear correlation between 
the frequency identified as v u , and V£ ow , Vh , and Vb 
(Fig. 18). This would not be expected if this fre- 
quency is really VhHz , as in the atoll sources VhHz 
remains approximately constant while the other 
frequencies change. Note also that with the ad- 
dition of the extreme island state results for Aql 
X-l, the atoll sources now also show these same 
correlations in the extreme island state (when v u is 
around 200 Hz). So, we identify the broad low fre- 
quency features previously reported by Wijnands 
& van dcr Klis (1998b) from the 1998 outburst as 
manifestations of L u . 

The broad Le ow component found in the ex- 
treme island states of the atoll sources and in the 
low luminosity bursters has now also been iden- 
tified in the accreting millisecond pulsars. It was 
suggested that this component could be identified 
with the lower kilohcrtz QPO based on extrapola- 
tions of frequency-frequency relations (see Psaltis 
et al. 1999a; Belloni et al. 2002; van Straaten et al. 
2002, 2003), however, it was also noted that it 
could be a different component (van Straaten et al. 
2003) . Here we observe that if we make the rela- 
tions of SAX J1808.4-3658 and XTE J0929-314 
coincide with those of the atoll sources by multi- 
plying v u with 1.454, the vi ow vs. v u relation of 
SAX J1808.4-3658 does indeed coincide with the 
atoll sources one, but the one point of SAX 
J1808.4-3658 is still offset from the atoll relation. 
We also have to multiply vg with a similar fac- 
tor of about 1.5 to make the Lg point of SAX 
J1808. 4-3658 fall on the atoll relation. If we apply 
the alternative shift (see §5.1) we have to multi- 
ply vg ow with about 0.35, but vt with a different 
factor of about 0.8 to make the relations coincide. 
Therefore it seems probable that Le ow and Lg are 
different components. 

For SAX J1808.4-3658 the characteristic fre- 
quency of LhHz seems to be slightly correlated 
with the characteristic frequencies of the other 
components (see Fig. 18). The LhHz points also 
seem to lie on an extension of the Lg ow relation 
in Figure 18. In the other atoll sources Li ow dis- 
appears and LhHz appears when the source goes 
from the extreme island to the island state. This 
happens at v u ~ 350 Hz in Figure 20. It might be 
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that in SAX J1808. 4-3658 Lg ow is present in the is- 
land state, but is not detected. It is not detected 
as a separate component when its characteristic 
frequency is in the range of VhHz, causing the cor- 
relation of LhHz with v u in Figure 18. However 
note that the scatter in the LhHz relation of the 
atoll sources is large and that a similar scatter in 
SAX J1808.4-3658 might accidentally cause the 
correlation of LhHz with v u . 

In the atoll sources 4U 0614+09, 4U 1608-52, 
and 4U 1728-34 the behavior of the band- limited 
noise in the lower left banana state is complex. 
First, the characteristic frequency of the band- 
limited noise increases from ~ 0.3 to ~ 15 Hz. At 
this stage the band-limited noise is broad. Then at 
~ 15 Hz this noise component appears to "trans- 
form" into a narrow QPO whose frequency in- 
creases up to ~ 50 Hz, while what appears to be 
another band-limited noise component appears at 
lower characteristic frequencies. In our previous 
papers (van Straaten et al. 2002, 2003) we call 
the band-limited noise component that becomes 
a QPO, as well as the QPO it becomes, L& and 
the "new" broad band-limited noise appearing at 
lower frequency L b2 . In SAX J1808.4-3658 the 
same phenomenon seems to occur in groups 2 and 
3. Based on this similar behavior the <~ 15 Hz 
Lorentzian can be identified as Lb and the broad 
component at ~ 10 Hz as Lb 2 . However, in Fig- 
ure 20 for v u > 600 Hz the L h relation of SAX 
J1808. 4-3658 follows the relation of the Lb 2 com- 
ponent of the atoll sources. The L(, 2 points of SAX 
J1808.4-3658 fall below that relation. This differ- 
ence between Lf, in SAX J1808. 4-3658 and in the 
other atoll sources can also be seen for z/& > 10 Hz 
in Figure 19. It seems that in this regime the fre- 
quency behavior of the band-limited noise in SAX 
J1808. 4-3658 deviates from that in the other atoll 
sources. 

In addition to the broad components Lb, Lh, 
Lg ow , the low luminosity bursters (Belloni et al. 

2002) and the atoll source 4U 1608-52 (van 
Straaten et al. 2003; Yoshida et al. 1993) show 
a narrow QPO that has f max between z/& and 
v h . The QPO in 4U 1608-52 is not the same as 
that in the low luminosity bursters (called Llf) 
but might be a sub-harmonic (van Straaten et al. 

2003) . Therefore we will call it L LF / 2 . Llf and 
Llf 72 are not only present in neutron star X-ray 
binaries but could also be identified in the black 



hole candidate GX 339-4 (Belloni et al. 2002; van 
Straaten et al. 2003). Also, a narrow QPO found 
in Cyg X-l by Pottschmidt et al. (2003) might be 
Llf 72 ( y an Straaten et al. 2003). Until now these 
QPOs were only found at relatively low frequencies 
(< 3 Hz). In this work we find a QPO with a v max 
between v\> and Vh at high frequencies (« 47 Hz) 
in SAX J1808.4-3658 (groups 2 and 3), and at low 
frequencies (0.7-1.0 Hz) in XTE J1814-338 (group 
29) and XTE J0929-314 (groups 23-25, previously 
reported by Markwardt et al. 2002). In Figure 22 
we plot our results of the SAX J1808.4-3658, XTE 
J1814-338, and XTE J0929-314 on Figure 13 of 
van Straaten et al. (2003) . Based on the power law 
relations of Figure 22 we can identify the QPOs 
in XTE J1814-338 and XTE J0929-314 as L LF/2 
and the QPOs in SAX J1808.4-3658 as Llf- 

5.3. Timing-Color Correlations and States 

The timing properties of the accreting millisec- 
ond pulsars SAX J1808.4-3658, XTE J1751-305, 
XTE J0929-314, and XTE J1814-338 are very 
similar to those of the atoll sources. Also based 
on its correlated timing and color behavior SAX 
J1808. 4-3658 can be classified as a classical atoll 
source (see §4.1). In the time span for which 
timing results could be obtained, XTE J0929-314 
and XTE J1814-338 showed little change in power 
spectra and colors. Therefore it is hard to clas- 
sify these sources as atoll sources although their 
correlated timing and color behavior is consistent 
with that of a typical atoll source (see §4.3 and 
§4.4). XTE J1751-305 is different from a classical 
atoll source with respect to the behavior of its hard 
color (see §3.3). In atoll sources in the island and 
extreme island states the characteristic frequen- 
cies and hard color are anti-correlated with each 
other, and decoupled from soft color and inten- 
sity, which are correlated. The hard color in XTE 
J 175 1-305 is instead strongly correlated with soft 
color and intensity, and is not anti-correlated with 
any of the characteristic frequencies. 

To explain the one-to-one relation between soft 
color and intensity in an extreme island state of 
the atoll source 4U 1608-52 van Straaten et al. 
(2003) used the following scenario: In the extreme 
island state of 4U 1608-52 the disk is far away from 
the neutron star and cannot be seen in the PCA 
energy spectrum; the soft color therefore only de- 
pends on the spectral properties of the neutron 
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star (Gierliriski & Done 2002b). As the accretion 
rate (intensity) increases, the neutron star surface 
becomes hotter and the soft color increases (van 
Straaten et al. 2003). This scenario can also be ap- 
plied to the horizontally extended extreme island 
state of SAX J1808.4-3658. In an energy spec- 
tral study of the 1998 outburst of SAX J1808.4- 
3658, Gierlihski et al. (2002) found that the soft 
component likely originates from the neutron star 
surface. 

In the horizontally extended island state branches 
that are observed in the atoll sources (including 
SAX J1808.4-3658), contrary to the horizontal 
branch of Z sources, the timing properties remain 
similar when the intensity changes. As noted by 
Olive et al. (2003), this aspect of the behavior can 
be explained by the scenario of van der Klis (2001) 
for the parallel track phenomenon in the intensity 
versus lower kilohertz QPO frequency diagram 
that occurs in the lower left banana of the atoll 
sources (see e.g. Mendez 2000). If we look at these 
island state branches in SAX J1808. 4-3658 during 
the 2002 outburst (groups 5-15) we observe two 
parallel tracks (see Figure 23), similar to those in 
the lower left banana of the atoll sources. This is 
a further indication that the same physical phe- 
nomenon is causing the source to drift between 
parallel tracks in the lower left banana of the atoll 
sources and along the island state branches as sug- 
gested by van Straaten et al. (2003). We observe 
one parallel track from group 5 to 8 and the next 
from group 9 to 13. Both parallel tracks last 4 
days; this is longer than the parallel tracks in the 
lower left bananas which usually last about a day. 

5.4. Ingredients for a Physical Model 

Currently, to our knowledge, there is no physi- 
cal model (see van Straaten et al. 2003, for a dis- 
cussion of models) that explains the mutual rela- 
tions between the characteristic frequencies of the 
power spectral features (e.g. Figures 18 and 19). 
Here we point out some clues that might help to 
construct such a model. As described above in 
§5.2, in the lower left banana of the atoll sources 
(among which we now also classify SAX J1808.4- 
3658) the behavior of the low frequency compo- 
nents is complex; several timing features, addi- 
tional to the ones observed in the island states, 
appear (see e.g. groups 2 and 3 of SAX J1808.4- 
3658). In Figure 20 this happens at v u « 600 



Hz. The Z sources also show a complex low fre- 
quency behavior (sec Figure 14 in van Straaten 
et al. 2003), also at high v u . If v u indeed represents 
the Keplerian frequency at the inner disk radius, 
a high v u means that then the inner disk is close 
to the neutron star surface, so the more complex 
behavior might be caused by relativistic or mag- 
netic effects that are stronger close to the neutron 
star. Note that the power spectra in the extreme 
island state of some sources also show extra timing 
features, but these are in the form of weak narrow 
QPOs that do not dominate the power spectra. 

Because of the complex behavior of the low- 
frequency noise above v u w 600 Hz in Figure 20 
we concentrate now on the more simple frequency 
relations below v u k> 600 Hz in Figure 20. The 
low-frequency relations below v u w 600 Hz can be 
fitted with power laws (see Table 7). For these 
fits we use the same method as in §5.1, which uses 
the FITEXY routine of Press Press & Teukolsky 
(1992). Note that although the x 2 /dof are rela- 
tively high, this is not due to systematic devia- 
tions from a power law but to scatter around it. 
The power law indexes of the fb, v^, and U( ow vs. 
v u relations fall in the range 1.8-2.7, progressively 
steeper towards lower frequency components. The 
errors on the power law indexes of the Llf and 
Llf/2 relations are large as these relations are only 
based on a few points, but also these power law in- 
dexes are consistent with this picture. 

Here we use v max as the characteristic fre- 
quency. Although f max has proven to be of great 
practical value and may represent a real underly- 
ing physical frequency (e.g., the inner frequency 
of a disk annulus, see Belloni et al. 2002), this 
is not assured. Two other frequencies, that have 
physical interpretations in terms of a damped har- 
monic oscillator, are vq, the centroid frequency 
of the Lorcntzian and ^damped (f° r a definition, 
see Titarchuk 2002). These frequencies can be 
calculated from our fit results ^ max and Q (see 
van Straaten et al. 2003). Note that using isq in- 
stead of ^ max does not give further insights in the 
shift between SAX J1808.4-3658, and the atoll 
sources (see §5.1). For the narrow Lorentzians 
the shift factor will be the same as z/ max is ap- 
proximately equal to vq, for the broad features vq 
is badly constrained and a shift cannot be deter- 
mined. For the v^.q vs. v u $ relation of the atoll 
sources 4U 0614+09, 4U 1608-52, and 4U 1728-34 
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van Straaten et al. (2003) found a quadratic re- 
lation. A quadratic relation was also found for 
the Vh& vs. u U fi relations of several Z sources 
(Stella & Vietri 1998; Psaltis et al. 1999b; Jonker 
et al. 1998; Homan et al. 2002). The relativistic 
precession model (Stella & Vietri 1998) predicts 
a quadratic relation between the Lcnsc-Thirring 
precession frequency and the Keplerian frequency 
of the inner disk; v LT = 4.4 x 10~ 8 I^M~ x v s v\ 
Hz, where -Z45 is the moment of inertia in units 
of 10 45 g cm 2 , M is the mass of the neutron star 
in Mq, and v s is the neutron star spin frequency. 
However, the required J45/M values were too large 
for hh to be the fundamental Lense-Thirring pre- 
cession frequency (Psaltis et al. 1999b). Also the 
Vh.o vs. v U fi relations coincided for 4U 1728-34 
and 4U 1608-52, which is not expected in this in- 
terpretation in view of their different spin frequen- 
cies (363 Hz and 620 Hz, respectively, assuming 
burst oscillations happen at the spin frequency). 
The relation of the frequency difference between 
the 410 Hz QPO and the ~ 401 Hz pulse spike of 
SAX J1808. 4-3658 was found to be consistent with 
a quadratic relation by Wijnands et al. (2003). We 
note, that this frequency difference is always about 
half of Vh,o (see Table 5). 

Although the relations between the character- 
istic frequencies of the low frequency components 
and v u are different between SAX J1808. 4-3658, 
XTE J0929-314 and the other sources (see §5.1), 
the mutual relations between the low frequency 
components are the same for all atoll sources (Fig. 
19). The Vh vs. i/j relation is even the same for 
black hole candidates and atoll sources (Wijnands 
& van der Klis 1998b; van Straaten et al. 2002). 
So, the Vh vs. Vb relation does not depend on 
spin, mass, magnetic held, or the presence of a 
surface. Note, that the frequency range of Vh and 
i>b is about a factor 10 lower for black hole candi- 
dates than for neutron star low mass X-ray bina- 
ries. Therefore and L/, might be properties of 
the accretion disk that scale inversely with mass 
(Wijnands & van der Klis 1999). 

The frequency correlations of SAX J1808.4- 
3658 and XTE J0929-314 are shifted compared to 
those of the other sources in a way that is most eas- 
ily described as a shift in upper and lower kilohertz 
QPO frequency by about a factor 1.5 (see §5.1). 
In most models v u reflects the Keplerian frequency 
at the inner edge of the accretion disk (e.g. Miller 



et al. (1998)). If this is indeed the correct interpre- 
tation, then for similar characteristic frequencies 
of the low frequency components, this inner radius 
of the disk can somehow get closer to the neu- 
tron star surface in the atoll sources, low luminos- 
ity bursters, and the accreting millisecond pulsars 
XTE J1751-305 and XTE J1814-338 than in the 
accreting millisecond pulsars SAX J1808. 4-3658 
and XTE J0929-314. It is unclear which physical 
parameter is the cause of this. In Table 8 we list 
several characteristic parameters for the accreting 
millisecond pulsars. There is no obvious relation 
between the shift and any of the parameters. In 
particular, there is no relation to spin frequency: 
SAX J1808.4-3658, and XTE J0929-314 which 
share one set of frequency relations have spin fre- 
quencies of respectively 401 and 185 Hz, and the 
other sources which share another set of relations 
have spin frequencies that range from 314 to 620 
Hz. Also the magnetic field strength B, which 
may be stronger for the accreting millisecond pul- 
sars than for the atoll sources (see e.g. Wijnands 
et al. 2003; Strohmayer et al. 2003), is not a good 
candidate as not all accreting millisecond pulsars 
show the shift. The fact that the shift occurs by 
a factor close to ~ 3/2 suggests a possible rela- 
tion with recently proposed parametric resonance 
models for kilohertz QPOs (e.g. Abramowicz et al. 
2003). In these models 2:3 frequency resonances 
between general relativistic orbital/epicyclic fre- 
quencies play a central role. However, similarly to 
the remarks about inner disk radius above, it is 
not clear what would set SAX J1808.4-3658 and 
XTE J0929-314 apart from the other sources. 
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Table 1 

Characteristic frequencies of the multi-Lorentzian fit for the accreting millisecond 



PULSARS AND LOW LUMINOSITY BURSTERS 



Group 
Number 


»m» (Hz) 


L t 2 
»m» (Hz) 


»m» (Hz) 


L IF L IF/2 
"max (Hz) (Hz) 


Wow 
»m» (Hz) 


'mo (Hz) 


»m» (Hz) 


L u 
"max (Hz) 


the accreting millisecond pulsars 


SAX .11808.4 :!(>r,H {the 2002 outburst) 


1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 


10.34±0.17 
15.66±0.28 
16.02±0.25 
7.204±0.097 
3.165±0.064 
3.083±0.058 
2.991±0.049 
2.264±0.056 
5.05±0.11 
3.94±0.11 
3.68±0.11 
3.51±0.10 
3.205±0.077 
3.93±0.13 
4.54±0.75 


8.7±1.0 
10.4±1.0 


48.4±1.2 
78.4±1.7 
74.9±1.2 
35.78±0.37 
18.21±0.19 
18.63±0.20 
18.17±0.19 
13.52±0.28 
25.41±0.28 
21.84±0.32 
20.30±0.34 
21.19±0.44 
19.44±0.30 
21.85±0.34 
30.0±2.0 


47. Oil. 1 
47.28±0.81 




196±17 
258±29 
329±42 
134.9±9.3 
120.6±9.1 
124±11 
149±24 
106±12 
118.6±8.7 
119±14 
110±12 
110±23 
148±25 
173±30 


503.6±5.3 


606.0±4.4 
691. 2 = 8.7 
685.1±5.6 
497.6±6.9 
352.1±6.5 
339.5±6.2 
328.5±9.2 
299.9±8.0 
396.4±4.2 
354.7±8.0 
349±12 
341±10 
337.5±7.6 
352±12 

374+ 48 
3 ' -29 


SAX J1808.4 3|j-,s ;tli. lulls outburst) 


16 
17 
18 
19 
20 
21 
22 


1.940±0.077 
0.407±0.011 
0.366±0.012 
0.3476±0.0093 
0.351±0.022 
0.987±0.035 
0.969±0.054 




12.44±0.37 
2.861±0.056 
2.612±0.064 
2.583±0.065 
2.77±0.18 
5.72±0.57 
2 




25.0±2.0 
21.5±1.9 
26.0±2.3 
27.8±8.8 
57.1±7.1 
2 


186 + 44 1 




269±17 1 
156.3±5.8 
151.2±6.8 
170.0±8.9 
169±22 
216±14 
2 


XTE .11751-305 


23 
24 


1.202±0.068 
0.897±0.067 
0.490±0.013 




6.10±0.16 
4.229±0.092 
2.942±0.034 




32.6±4.7 
23.0±2.6 
20.1±1.4 






278±13 
260±56 
234±12 


XTE .10929-314 


26 
27 
28 


0.250±0.014 
0.463±0.048 
0.42±0.10 




2.01±0. 17 
3.08±0.14 
2.65±0.40 


0.742±0.013 
1.072±0.022 
0.846±0.046 


31.4±4.7 

_3 
_3 






149±17 

_3 
_3 


XTE .11814-338 


29 
30 


0.381±0.015 
0.500±0.026 




2.173±0.050 
3.048±0.099 


0.910±0.031 


21.0±6.2 

49+ 31 
"-11 






191±30 
220±f 9 


the low luminosity bursters 


IE 1724-3045 


C 
D 
E 
F 
G 
H 
I 
J 
K 
L 


0.1614±0.0034 
0.555±0.035 
0.311±0.026 
0.1338±0.0074 
0.1392±0.0075 
0.1276±0.0048 
0.1362±0.0066 
0.1647±0.0073 
0.207±0.014 
0.175±0.011 




1.0206±0.0064 
3.421±0.100 
1.936±0.048 
0.822±0.018 
0.905±0.021 
0.746±0.015 
0.879±0.025 
1.024±0.017 
1.391±0.058 
1.155±0.046 


0.6246±0.0049 
0.6854±0.0082 
0.5137±0.0048 

°- 888 _0.007 
0.7676±0.0091 
1.015±0.015 
0.877±0.021 


11.12±0.21 
17.4±3.5 
16.8±1.3 
9.65±0.25 
9.75±0.37 
8.54±0.21 
9.67±0.45 
11.27±0.30 
12.43±0.80 
12.62±0.61 






156.0±6.3 
230±26 
190±17 
146±10 
138±10 

133.7±7.6 
169±23 

167.7±9.4 
157±21 
173±17 


SLX 1735-269 


M 


0.1642±0.0090 




1.155±0.033 




13.5±1.6 






176±26 


GS 1826-24 


N 
O 
P 
Q 

R 


0.264±0.011 
0.2208±0.0081 
0.437±0.032 
0.805±0.043 
0.568±0.033 




1.787±0.040 
1.525±0.031 
2.609±0.084 
4.86±0.16 
3.79±0.31 


1.128±0.022 
0.975±0.020 
1.814±0.039 

2 3Q + 0- 11 
2 ' JJ -0.07 


15.51±0.67 
14.42±0.49 
21.2±1.2 
24.8±2.7 
20.87±0.98 






178±24 
167±16 
187 = 31 
292 ±50 
240±35 


1 the power 


of this Lorentzian 


could only be t 


neasured to an a 


ccuracy of better than 50% (and not 


33%, see §3.2). 








2 due to bad 


statistics the po' 


wer above 4 Hz 


was fitted with e 


inly one broad Lorentzian. L^, ^£ OVJ 


, and L u could 


not be identified (see §3.2). 





the high frequency part of the power spectra is fitted with one broad Lorentzian that probably fits both L^ ow and L u (sec §3.4). 



NOTE. — isted are the characteristic frequencies (= fmax) of the different Lorcntzians used to fit the power spectra of the accreting millisecond pulsars (described 
in §3) and the low luminosity bursters (see §5.1). The groups of the low luminosity bursters are the same as in Belloni et al. (2002). The quoted errors in i/ max 
use Ax 2 = 1-0 (corresponding to a 68% confidence level). 
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Table 2 

q values of the multi-lorentzian fit for the accreting millisecond pulsars and low 

luminosity bursters 



Group 


Lb 


L b2 








Lf 




Number 


Q 


Q 


Q 


Q Q Q 


Q 


Q 


Q 


the accreting millisecond pulsars 


SAX .11808.4- acf,K ;tlu' 2002 outburst) 


1 


0.322 + 0.013 




1.03±0.11 




0.60±0.26 




5.9±1.0 


2 


1.36±0.20 


0.337±0.034 


2.47±0.57 


2.25±0.37 


0.83 + 0.38 




5.17±0.95 


3 


, „,+0.28 


0.272 + 0.023 


2.66±0.52 


2.38±0.30 


0.35±0.18 


14.26 (fixed) 


5.60±0.77 


4 


0.312±0.012 




1.195 + 0.061 




0.43±0.12 




2.75±0.32 


5 


0.238±0.015 




0.813±0.042 




0.281±0.094 




2.12±0.30 


6 


0.257 + 0.014 




0.807±0.042 




0.28±0.11 




2.22 + 0.35 


7 


0.286±0.014 




0.914±0.050 




0.03±0.17 




3 1 + 1 - 6 

d -0.8 


8 


0.211±0.017 




0.642 + 0.065 




0.18±0.10 




1.70±0.42 


9 


0.255±0.016 




1.136±0.073 




0.32±0.13 




4.13±0.55 


10 


0.234 + 0.019 




0.890 + 0.065 




0.30±0.17 




2.63±0.54 


11 


0.245±0.021 




0.949±0.088 




0.34±0.19 




, ,,+0.68 
2 - -0.45 


12 


0.217±0.021 




0.875±0.087 




0.35±0.30 




2 07+ ' 78 
• iu '-0.52 


13 


0.237±0.017 




0.891 + 0.067 




0.22±0.15 




„ -.4-0.99 
' -0.69 


14 


0.208±0.022 




1.085±0.090 




(fixed) 




3 4+ 1 ' 3 
3 -0.9 


15 


0.13±0.10 




0.85±0.26 








1.47+0.77 



SAX Jisns. I :',ii.-.8 (the 19 



16 
17 
18 
19 
20 
21 
22 


0.342 + 0.035 0.489 + 0.074 
0.186 + 0.022 0.413 + 0.045 
0.167 + 0.027 0.450 + 0.061 
0.222 + 0.021 0.240 + 0.050 
0.309 + 0.060 0.18 + 0.13 
0.186 + 0.020 (fixed) 
0.020 + 0.045 - 2 




(fixed) 
(fixed) 
(fixed) 
(fixed) 
0.28 + 0.16 
_2 


(fixed) 1 


1 6+ 1 -' i 1 
0.362 + 0.064 
0.332 + 0.075 
0.278 + 0.077 
0.26 + 0.26 
0.54+0.13 
2 


XTE .11751-305 


23 
24 
25 


0.167+0.033 0.70 + 0.10 
0.106 + 0.043 0.787 + 0.091 
0.057 + 0.018 0.724 + 0.038 




13 + u - 2:i 
u l ' -0.93 

0.39+0.20 

(fixed) 




1.01 + 0.28 
0.08+"; 22 
0.70 + 0.12 


XTE .10929-314 


26 
27 
28 


0.76+0.15 (fixed) 
0.382 + 0.058 0.367 + 0.097 
0.49 + jJ 2 ° 0.11 + 0.21 


2.67 + 0.54 
, ,,. + 0.77 
2 - 30 -0.51 

2U -0.7 


0.47 + 0.20 
3 

_ 3 




0.43 + 0.21 

_3 

_3 


XTE .11814-338 


29 
30 


0.251 + 0.026 0.581 + 0.071 
0.335 + 0.037 0.290 + 0.065 


3 5 + 1 - 5 


(fixed) 

13+ - 20 
ui -0.31 




09+ U iU 

-0.33 

41 + 33 
u41 -0.15 


the low luminosity bursters 


IE 1724-3045 


C 
D 
E 
F 


0.122±0.019 0.570+0.016 
0.139 + 0.044 0.653 + 0.093 
0.117 + 0.059 0.600 + 0.083 
0.135 + 0.043 0.470 + 0.037 


7 , + 2.2 


(fixed) 
(fixed) 
0.10 + 0.13 
0.081 + 0.042 




0.268±0.051 
0.38 + 0.15 
0.72 + 0.20 

0.208 + 0.085 


G 
H 
I 
J 

K 


0.148 + 0.044 0.512 + 0.044 
(fixed) 0.535 + 0.032 
0.210 + 0.046 0.479 + 0.057 
0.083 + 0.036 0.553 + 0.031 
0.122 + 0.049 0.505 + 0.064 


r . + 8.9 
— 2 2 
6.21 + 0.82 
7 +17 
-3 
7.1 + 1.3 

4 4 + 1 ' 6 
*■ — 1.1 


0.036 + 0.068 
(fixed) 

0.033 + 0.085 
(fixed) 
0.07 + 0.12 




0.30 + 0.12 
0.229 + 0.078 

11 + - 15 
-0.25 
0.375+0.078 
0.25 + 0.16 


L 


0.111 + 0.045 0.491 + 0.050 


3 2 + 1.6 
d - -0.9 


(fixed) 




0.30 + 0.13 


SLX 1735-269 


M 


0.040±0.042 0.457±0.061 




(fixed) 




0.33+0.16 


GS 1826-24 


N 
O 


0.050+0.029 0.562±0.050 
0.139 + 0.028 0.486+0.037 


7.1+2.9 

6 2+ 2 ' 5 
D -1.5 


0.082±0.075 
(fixed) 




0.38+0.23 
0.47 + 0.14 


P 
Q 

R 


0.049 + 0.045 0.69 + 0.10 
0.133 + 0.040 0.76 + 0.12 
0.017+0.043 0.684+0.095 


5 fi+ 2 - 6 
2.9+1.1 


0.05 + 0.10 
0.29 + 0.20 
0.43 + 0.11 




1 4 + 1 - 3 
-0.6 
0.82 + 0.44 
0.46 + 0.20 


1 thc 


power of this Lorentzian could only be measured to an 


accuracy of better than 50% (and not 33%, sc 


M §3.2). 




2 due 


to bad statistics the power above 4 Hz was fitted with 


only one broad Lorentzian 




x could not be identified (sec 


§3.2). 


3 the 


high frequency part of the power spectra is fitted with 


one broad Lorentzian that 


probably fits both 


L fora » nd L « < sco § 3 ' 4 )- 





NOTE. — isted are the Q values {= ^q/2A) of the different Lorentzians used to fit the power spectra of the accreting millisecond pulsars (described in §3) 
and the low luminosity bursters (see §5.1). The groups of the low luminosity bursters are the same as in Bclloni ct al. (2002). The quoted errors in Q use 
Ax 2 = 1.0 (corresponding to a 68% confidence level). 
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Table 3 

Integrated fractional rms of the multi-Lorentzian fit for the accreting millisecond 



PULSARS AND LOW LUMINOSITY BURSTERS 



Group 


H 


L b2 




L LF L LF/2 




^hHz 




L, 




rms (%) 


rms (%) 


rms (%) 


rms (%) rms (%) 


rms (%) 


rms (%) 


rms (%) 


rms (%) 


the accreting millisecond pulsars 


SAX .J180S.1 :M).-,R ; tli. 2002 outburst) 


1 


17.05±0.14 




11.78±0.65 






10.2 + 1.1 




8.92 + 0.50 


2 


8.68±0.91 


10.13±0.85 


6.96±0.75 


6.87±0.57 




7.98+0.98 




8.32 + 0.57 


3 


8.7±1.2 


12.55±0.96 


6.79±0.66 


7.17±0.49 




11.66+0.95 


2.94 + 0.48 


8.78 + 0.58 


4 


16.61±0.11 




13.94±0.38 






12.43+0.72 




10.87+0.49 


5 


15.40±0.15 




15.71±0.35 






14.99+0.78 




10.84 + 0.74 


6 


15.11±0.14 




15.73±0.35 






14.71 + 0.87 




10.76 + 0.82 


7 


15.23±0.12 




15.25±0.40 






17.3 + 1.3 




8.8+1.5 


8 


16.23±0.20 




14.21 + 0.59 






17.0+1.1 




11.4+1.3 


9 


16.03±0.16 




14.58±0.48 






14.39+0.87 




10.29 + 0.57 


10 


15.30±0.21 




16.11 + 0.56 






14.1 + 1.2 




10.67+0.99 


11 


15.56±0.22 




15.58±0.67 






14.5 + 1.4 




10.9+1.2 


12 


15.67±0.22 




16.03 + 0.84 






13.9 + 2.2 




11.9+1.6 


13 


16.37±0.18 




15.60±0.51 






16.5 + 1.5 




10.7+1.5 


14 


16.14±0.23 




14.28±0.48 






19.22 + 0.93 




9.0+1.4 


15 


15.4±1.1 




17.1±1.5 






< 16. 1 




19.5 + 2.8 


SAX J1808.4 3658 (the 1998 outburst) 


16 
17 


15.32 + 0.38 
19.76±0.25 




16.09±0.61 
15.85±0.49 




18.14+0.46 


18.7+-J ^ 




10 i + 9 - 5 2 

-2.6 
19.85 + 0.68 


18 


20.85±0.31 




15.34+0.64 




19.04+0.51 






20.58 + 0.79 


19 


20.26±0.31 




17.33±0.60 




19.04+0.55 






19.65 + 0.85 


20 


18.83±0.82 




18. 0±1. 5 




18.1 + 2.0 






20.8 + 2.6 


21 


16.33±0.39 




16.53±0.35 




15.8+2.0 






17.3+1.6 


22 


25.08±0.48 




_3 




_3 






_3 


XTE J1751-305 


23 


10.70 + 0.29 




12 + "- 8 




11.2 + 2.5 






12.2+1.1 


24 


11.89±0.39 




12.37±0.65 




10.9 + 1.5 






15.3+1.5 


25 


12.53±0.13 




12.05±0.26 




14.57 + 0.23 






13.49 + 0.57 


XTE J0929-314 


26 


5.36±0.44 




15.99±0.25 


5.00±0.49 


13.3 + 2.1 






18.8 + 2.0 


27 


8.13±0.52 




12.66±0.61 


5.23±0.73 


_4 






_4 


28 


7.1±1.6 




12.6±1.2 


ou -1.2 


_4 






_4 


XTE J1814-338 


29 


11.82±0.28 




13.92±0.65 


2.89±0.58 


15.4+1.9 






22.0 + 2.1 


30 


9.83±0.46 




15.81 + 0.705 




ie..+«;5 








the low luminosity bursters 


IE 1724-3045 


C 


14.16±0.14 




16.06±0.14 




18.201 + 0.098 






13.87+0.30 


D 


12.00 + 0.33 




12.20±0.74 




12.24 + 0.52 






14.37+0.93 


E 


11.87±0.42 




13.57 + 0.69 




14.14 + 0.79 






11.73 + 0.92 


F 


14.37±0.39 




16.46±0.46 


3.25±0.31 


18.50 + 0.40 






14.75 + 0.56 


G 
H 


14.32±0.36 
16.07±0.24 




16.27±0.49 
15.15 + 0.35 


, r ,+0.66 
-0.32 
4.21 + 0.31 


18.20 + 0.57 
19.33+0.13 






14.13 + 0.76 
13.70 + 0.41 


I 
J 


14.12±0.38 
14.71±0.29 




16.00±0.55 
15.42±0.34 


„ ,„+0.74 
3.32±0.34 


18 47 +°'65 
-0.93 
18.39 + 0.13 






14.3+1.0 
13.65 + 0.45 


K 


14.29±0.44 




15.36±0.74 


4 47 +0.88 
-0.61 


16.71+0.99 






14.0+1.1 


L 


14.15±0.39 




15.38±0.59 


47 +1.3 
-0.8 


18.03 + 0.24 






13.39 + 0.71 


SLX 1735-269 


M 


9.95±0.24 




9.96±0.38 




12.37 + 0.41 






12.78 + 0.94 


GS 1826-24 


N 


12.66±0.21 




11.90±0.45 


1.92 + 0.52 


15.28 + 0.61 






9.01 + 0.99 


O 


14.28±0.24 




14.20±0.36 


2.38 + 0.39 


18.60 + 0.16 






10.27 + 0.63 


P 


14.14±0.41 




12.73+0.88 


3.15 + 0.57 


17.24+0.66 






7.5+1.4 


Q 


13.61±0.31 




13.36±1.00 




13.2+\\ 






10.0+1.5 


R 


14.25±0.34 




13.02±0.84 


* -0.7 


13.75 + 0.78 






12.6+1.1 



^-95 % confidence upper limits (sec §3.1). 

2 the power of this Lorcntzian could only be measured to an accuracy of better than 50% (and not 33%, see §3.2). 

3 duc to bad statistics the power above 4 Hz was fitted with only one broad Lorentzian. , ^>£ ow , and L u could not be identified (sec §3.2). 
4 the high frequency part of the power spectra is fitted with one broad Lorentzian that probably fits both L'£ ow and L u (see §3.4). 



NOTE. — isted are the values of the integrated fractional rms (over the full PCA energy band) of the different Lorcntzians used to fit the power 
Spectra of the accreting millisecond pulsars (described in §3) and the low luminosity bursters {see §5.1). The groups of the low luminosity bursters 
are the same as in Bclloni ct al. (2002). The quoted errors in the rms use A% 2 = 1.0 (corresponding to a 68% confidence level). 
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Table 4 




RXTE OBSERVATION IDS OF THE FIT GROUPS 


Group 


RXTE observation IDs 




n ri4 no A ni non 


2 


a m no a m A ni ru 


3 


A-01-01 


4 


a ni 090 a m 09 


5 


A 09 01 A 09 000 A 09 00 


g 


A 09 09 A 09 03 A 09 04 


7 


A-02-05 A-02-06 A-02-07 


8 


A-02-10 A-02-08 

A U^i 1 V 7 . A U U 


9 


A 09 90 A 09 0Q A 09 1 Q A 09 93 


10 


A-02-21 A-02-11 A-02-18 A-02-12 A-02-22 


H 


A-02-13 A-02-15 A-02-14 A-02-16 


12 


A 09 1 7 A 03 000 A 03 00 A 03 04 A 03 ft^ 

A Ui 1 I - A UO UUU, A UO UU • A UO Ut: • A UO UJ 


13 


A 03 Ofi A 03 07 A 03 010 A 03 01 A 03 08 A 03 0Q A 03 1 A 03 1 1 


14 


A-03-1 9 A-03-1 3 A-03-09 A-03-1 4 A-03-03 C-0V00 R-01 -04 


15 


R-01 -000 R-01 -00 

1J U 1 U U U i U U± UU 


16 


ri-0i-03S 

17 U 1 UJ J 


17 


11-03-00 D-04-00 


18 


D-05-00 

U U'7 UU 


19 


D-Ofi-OI D-06-000 D-06-00 

J— ' UU Ul) L/ UU UUU , 17 UU UU 


20 


D-07-00 


21 


r)_0S_00 D-OQ-01 D-OQ-09 D-0Q-03 D-OQ-04 

17 UU UU, 17 UJ Ulj 17 UJ Uij J7 UJ UJ) 17 UJ U^ 


22 


d oq oo n io 09 n io 01 n 10 00 

.17 UJ UU) 17 J.U UZi, 17 J.U U-L, .17 -1U UU 


23 


F 01 00 F 01 00 F 09 00 

J_j U -1 UU) -1 U -1 UU) 1_J Ui UU 


24 


F-02-00 


25 


F-03-01, F-03-00, F-04-00 F-05-03, F-05-02, F-05-01, F-05-04, F-05-000, F-05-00, F-06-00, 




F-06-01, F-07-00, F-07-01, F-08-000, F-08-00, F-09-000 


26 


G-03-01, G-03-00, G-04-00, G-04-01, G-05-01, G-07-01 


27 


G-02-00, G-02-01, G-06-01, G-06-00, G-05-00, G-07-00 


28 


G-08-00, G-08-01, G-08-02, G-08-03, G-08-04, G-08-05, G-08-06, G-09-00, G-09-01, G-09-02, 




G-10-00, G-10-01, G-10-02, G-10-03, G-10-04, G-11-00, G-12-00, G-12-01, G-12-02 


29 


H-01-01, H-01-02, H-01-03, H-01-04, H-01-05, H-01-06, H-02-04, H-02-05, H-03-12, H-03-13, 




H-03-000, H-03-00, H-03-01, H-03-03, H-03-04, H-03-05, H-03-02, H-03-06, H-03-07 


30 


H-01-07, H-01-08, H-01-09, H-02-01, H-02-00, H-02-03, H-02-02, H-02-06, H-02-07, H-02-09, 




H-02-08, H-03-08, H-03-09, H-03-10, H-03-11, H-04-07, H-04-00, H-04-01, H-04-02, H-04-08, 




H-04-03, H-04-04, H-04-05, H-04-09, H-04-06, H-05-00, H-05-01 



NOTE. — Listed arc RXTE observation IDs of each fit group. The letters stand for part of the observation ID; 
A = 70080-01, B = 70080-02, C = 70080-03, D = 30411-01. E = 70134-03, F = 70131-01, G = 70096-03, and H 
= 80418-01. 

Table 5 

Characteristic frequencies of the 410 Hz QPO 



Groups 


410 Hz QPO 
tWx (Hz) 


^410 - ^pulsc 
(Hz) 


tWx (Hz) 


vo (Hz) 


5-6 
9 

10-11 


409.2t ( i'o 
413.41±0.30 
409.3±1.7 


8.3tr 

12.47±0.30 
8.36±1.7 


18.67±0.15 
25.55±0.29 
20.95±0.24 


16.81±0.13 
23.41±0.26 
18.54±0.21 



Note. — isted are the characteristic frequency (i/ max ) of the 410 Hz 
QPO, the frequency difference between the ~ 401 Hz pulse spike (^ pu i se ) 
and the 410 Hz QPO (^410 — '-'pulse)- Also listed are the i/ maI and vq of 
L h . U() is the centroid frequency of the Lorcntzian and can be obtained 
from our fit parameters /-'max and Q as uq = 2K max Q/ v^Q 2 + 1- For the 
410 Hz QPO i/max = fo as Q is between 50 and 100. 
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Table 6 






Shift factors 






v u shift factor Vi shift factor 
(u u X factor) (vi X factor) 


X 2 /dof 1 


SAX J1808.4-3658 


vt 


1.420±0.013 0.355±0.013 
1.481±0.013 0.375±0.011 
1.653±0.073 0.334±0.061 
1.359±0.021 
0.804±0.015 


35.4/31 
54.5/31 
7.8/9 
183.8/15 
183.8/15 


XTE J0929-314 


fb 

v tow 


1.40±0.16 0.35±0.13 
1.41±0.16 0.42±0.13 
2.03±0.30 0.283±0.090 


15.8/15 
30.1/15 
5.5/5 


XTE J1751-305 




1.188±0.045 0.585±0.073 
1.112±0.042 0.763±0.076 
1.036±0.058 0.93±0.10 


18.6/17 
32.9/17 
6.3/7 


XTE J1814-338 


Vb 
Vh 


1.23±0.14 0.53±0.19 
1.13±0.13 0.73±0.22 
1.42±0.24 0.52±0.16 


15.8/16 
30.1/16 
6.5/6 


IE 1724-3045 


Vb 
Vh 


1.167±0.043 0.622±0.076 
1.037±0.040 0.915±0.091 
1.051±0.077 0.92±0.12 


22.2/24 
36.7/24 
9.7/14 


SLX 1735-269 


»b 
Vh 
Vlow 


1.04±0.15 0.90±0.40 
0.96±0.14 1.11±0.41 
1.07±0.18 0.88±0.27 


15.8/15 
30.1/15 
5.5/5 


GS 1826-24 


Vb 
Vh 
v low 


1.195±0.071 0.58±0.11 
1.134±0.067 0.73±0.11 
1.117±0.075 0.83±0.11 


17.4/19 
31.1/19 
9.1/9 



1 the quoted x 2 /dof are independent of whether the 
shift was in Vi or v u . 

2 after v u was multiplied with 1.454 (see §5.1). 

Note. — Shift factors between the frequency relations 
of the listed sources and the atoll sources (see §5.1). The 
errors in the shift factor use A% 2 = 1.0. 

Table 7 



Power law fits to the frequency relations 



Component 


Index 


Normalization 


X 2 /d°f 


u 


2.724±0.038 


1.49(±0.34) X 10~ 7 


152/53 




2.439±0.030 


4.62(±0.85) x 10~ 6 


89/53 




1.81±0.11 


1.15(±0.66) x 10~ 3 


34/31 




3.00±0.63 


2.0(±6.9) x 10~ 7 


6/9 


Llf/2 


3.13±0.84 


4(±22) x 10~ 8 


1.7/1 



Note. — isted are the index, the normalization and the x 2 /dof 
of power law fits to the i/ max versus v u relations in Figure 20. 
Only the data below about v u fa 600 Hz is included in the fit. 
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Table 8 

Characteristic parameters for the accreting millisecond pulsars 



Source 


^ 'spin 


P orb 


Bursts 


Shift factor 


References 1 


SAX J1808. 4-3658 


401 Hz 


2 hr 


Y 


~ 1.5 


a, b, c, d 


XTE J0929-314 


185 Hz 


44 min 


N 


~ 1.5 


e 


XTE J1751-305 


435 Hz 


42 min 


N 2 


~ 1 


f 


XTE J1814-338 


314 Hz 


4.3 hr 


Y 


~ 1 


g, h, i 



x a: Wijnands & van der Klis (1998a); b: Chakrabarty & Morgan (1998); c: 
in't Zand ct al. (2001); d: Chakrabarty et al. (2003); e: Galloway ct al. (2002); 
f: Markwardt et al. (2002); g: Markwardt & Swank (2003); h: Markwardt et al. 
(2003); i: Strohmayer et al. (2003) 

2 the RXTE data of XTE J1751-305 shows 1 type I X-ray bursts which originated 
from from GRS 1747-312, see in't Zand ct al. (2003) 

Note. — isted are the source name, the spin frequency, v sp i n , the orbital period, 
Porf>> whether or not the source has shown type I X-ray bursts, and the shift factor 
that has to be applied to v u to make the relations coincide with those of the atoll 
sources (see §5.1). 
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0.1 10 1000 
Frequency (Hz) 

Fig. 1. — The multi-Lorentzian components of the Ijgv mass X-ray binary scheme of Bclloni ct al. (2002) 
and van Straaten et al. (2003). The states in which the components occur are indicated; lower left banana 
(LLB), island state (IS), and extreme island state (EIS). Note that Llf and are not shown. They 

generally occur in the extreme island state as a narrow QPO with a frequency between that of Lb and L/j. 
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Fig. 2. — Six representative power spectra and fit functions in the power spectral density times frequency 
representation for the 2002 outburst of SAX J1808. 4-3658. The different lines in the first 5 panels mark 
the individual Lorentzian components of the fit that are also indicated in the plots. The dashed lines in 
the sixth panel mark the Gaussian components of the fit, the dotted lines mark the Lorentzian components. 
The arrows indicate the 40 f Hz pulsar spike that was excluded during the fit. The group numbers are also 
indicated in the plots. 
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Fig. 3. — Hard color (top panel), soft color (second panel from the top), intensity (third panel from the 
top), and the characteristic frequencies of two of the power spectral components (bottom panel) plotted 
versus time for the 2002 outburst of SAX J1808. 4-3658. In the bottom panel the triangles mark Vh and the 
stars the centroid frequency of the main Gaussian peak used to fit the 1 Hz complex (see §3.1). For clarity 
some of the group numbers are indicated. Errors on time in the bottom plot indicate the addition of several 
observations to improve the statistics (see §2.2). 
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Fig. 4. — The combined power spectrum and fit function of groups 5-6 and 9-11 in the power spectral 
density times frequency representation. The 410 Hz QPO and the ~ 401 Hz pulse spike, that was excluded 
during the fit, arc indicated. 
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Fig. 5. — Three representative power spectra and fit functions in the power spectral density times frequency 
representation for the 1998 outburst of SAX J1808. 4-3658. The different lines mark the individual Lorcntzian 
components of the fit that are also indicated in the plots. The arrows indicate the 401 Hz pulsar spike that 
was excluded during the fit. The group numbers are also indicated in the plots. 
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Fig. 6. — Hard color (top panel), soft color (second panel from the top), intensity (third panel from the 
top), and v\, (bottom panel) plotted versus time for the 1998 outburst of SAX J1808. 4-3658. For clarity 
some of the group numbers are indicated. Errors on time in the bottom plot indicate the addition of several 
observations to improve the statistics (see §2.2). 
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Fig. 7. — Three representative power spectra and fit functions in the power spectral density times frequency 
representation of XTE J 175 1-305. The different lines mark the individual Lorentzian components of the fit 
that are also indicated in the plots. The arrows indicate the 435 Hz pulsar spike that was excluded during 
the fit. The group numbers are also indicated in the plots. 
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Fig. 8. — Hard color (top panel), soft color (second panel from the top), intensity (third panel from the top), 
and Vb (bottom panel) plotted versus time for the 2002 outburst of XTE J1751-305. Group numbers are 
indicated. Errors on time in the bottom plot indicate the addition of several observations to improve the 
statistics (see §2.2). 
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Fig. 9. — Three representative power spectra and fit functions in the power spectral density times frequency 
representation of XTE J0929 314. The different lines mark the individual Lorentzian components of the fit 
that are also indicated in the plots. The arrows indicate the 185 Hz pulsar spike that was excluded during 
the fit. The group numbers are also indicated in the plots. 
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Fig. 10. — Hard color (top panel), soft color (second panel from the top), intensity (third panel from the 
top), and Vb (bottom panel) plotted versus time for the 2002 outburst of XTE J0929-314. Errors on time in 
the bottom plot indicate the addition of several observations to improve the statistics (see §2.2). 
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Fig. 11. — Two representative power spectra and fit functions in the power spectral density times frequency 
representation of XTE J1814-338. The different lines mark the individual Lorentzian components of the fit 
that are also indicated in the plots. The arrows indicate the 314 Hz pulsar spike that was excluded during 
the fit. The group numbers are also indicated in the plots. 
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Fig. 12. — Hard color (top panel), soft color (second panel from the top), intensity (third panel from the 
top), and Vb (bottom panel) plotted versus time for the 2003 outburst of XTE J1814-338. Errors on time in 
the bottom plot indicate the addition of several observations to improve the statistics (see §2.2). 
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Fig. 13. — Color-color diagram (left panel) and soft color vs. intensity diagram (right panel) of SAX J1808.4- 
3658. The 1998 outburst is marked by open triangles, the 2002 outburst by filled squares. The line connects 
observations that are adjacent in time. The group numbers of the power spectra are indicated. Colors are 
in units of Crab (see §2). 
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Fig. 14.— Hard color plotted versus v h for SAX J1808.4-3658, XTE J1751-305, XTE J0929-314, and XTE 
J1814-338. The sources are indicated in the plot. 
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Fig. 15. — Color-color diagram (left panel) and soft color vs. intensity diagram (right panel) of XTE J1751- 
305. The black point mark the 2002 outburst and the reflare, the grey points mark the background points 
that are due to nearby sources and Galactic diffuse emission (see §3.3) The line connects observations that 
are adjacent in time. The group numbers of the power spectra are indicated. Colors are in units of Crab 
(see §2). 
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Fig. 16. — Color-color diagram (left panel) and soft color vs. intensity diagram (right panel) of XTE J0929- 
314. The filled squares mark group 26, the open diamonds group 27, and the filled stars group 28. The line 
connects observations that are adjacent in time. Colors are in units of Crab (see §2). 



38 




Soft color (Crab) Intensity (Crab) 

Fig. 17. — Color-color diagram (left panel) and soft color vs. intensity diagram (right panel) of XTE J1814- 
338. The filled squares mark group 29, the open diamonds group 30. The line connects observations that 
are adjacent in time. The insert in both panels is a blow up of that part of the diagrams for which we also 
obtained timing results. Colors are in units of Crab (see §2). 
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Fig. 18. — Correlations between the characteristic frequencies (= ^ max ) of the power spectral components 
plotted versus v u . The colored points mark the accreting millisecond pulsars SAX J1808. 4-3658 (red), XTE 
J1751-305 (blue), XTE J0929-314 (yellow), and XTE J1814-338 (green). The black point mark the atoll 
sources 4U 0614+09, 4U 1608-52, 4U 1728-34, and Aql X-l. The grey point mark the low luminosity bursters 
IE 1724-3045, GS 1826-24, and SLX 1735-269. The solid dots mark L b2 , the triangles L b , the crosses L h , 
the squares h^Hz, the stars Li and the diamonds Lg ow . 
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Fig. 19. — The characteristic frequencies vt ow (top panel) and Vh (bottom panel) plotted versus Vb- The 
filled grey squares mark the atoll sources 4U 0614+09, 4U 1608-52, 4U 1728-34, and Aql X-l, the open black 
diamonds mark the low luminosity bursters IE 1724-3045, GS 1826-24, and SLX 1735-269, and the filled 
black stars mark the accreting millisecond pulsars SAX J1808.4-3658, XTE J1751 305, XTE J0929-314, 
and XTE J1814-338. 
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Fig. 20.— As Figure 18, but with v u and v t of SAX J1808.4-3658 and XTE J0929-314 multiplied by a 
factor 1.454. 
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Fig. 21. — Q versus z/ ma x for L/, and L u . The black points mark SAX J1808. 4-3658, the gray points the atoll 
sources. 
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Fig. 22. — Characteristic frequency of versus that of the narrow low-frequency QPOs that have a ^ max 
between v\, and Vh- The sources are indicated in the plot. Note that we use the results of this paper for the 
low luminosity bursters IE 1724-3045 and GS 1826-24, which for the low frequencies are in agreement with 
the results of Belloni et al. (2002) that were used in Figure 13 of van Straaten et al. (2003). The dashed 
line indicates a power law fit to the vlf vs. Uh relation of the low-luminosity bursters IE 1724-3045 and GS 
1826-24, and the BHC GX 339-4. The dash-dotted line is a power law with a normalization half of that of 
the dashed line. 
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Fig. 23. — Parallel lines during an island state of SAX J1808. 4-3658. The top panel shows the v u vs. intensity 
diagram and the bottom panel the hard color vs. intensity diagram. The group numbers are indicated. The 
points that span up the parallel tracks are connected by a line. 
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